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ABSTRACT 

It is argued and demonstrated by particle-in-cell simulations that the syn- 
chrotron maser instability could develop at the front of a relativistic, magnetized 
shock. The instability generates strong low-frequency electromagnetic waves 
propagating both upstream and downstream of the shock. Upstream of the 
shock, these waves make electrons lag behind ions so that a longitudinal elec- 
tric field arises and the electrons are accelerated up to the ion kinetic energy. 
Then thermalization at the shock front results in a plasma with equal temper- 
atures of electrons and ions. Downstream of the shock, the amplitude of the 
maser-generated wave may exceed the strength of the shock-compressed back- 
ground magnetic field. In this case the shock-accelerated particles radiate via 
nonlinear Compton scattering rather than via a synchrotron mechanism. The 
spectrum of the radiation differs, in the low-frequency band, from that of the 
synchrotron radiation, providing possible observational tests of the model. 

Subject headings: instabilities-magnetic fields -masers-radiation mechanisms: 
non-thermal -shock waves 



1. Introduction. 



Relativistic shocks are supposed to play an important role in various astrophysical ob- 
jects, e.g., in pulsar wind nebulae, active galactic nuclei, gamma-ray bursts. In all these cases 
shocks are collisionless therefore some plasma instabilities are assumed to provide dissipation 
necessary for the final thermalization of the flow. An important point is that in relativistic 
case, plasma instabilities are generally electromagnetic and therefore strong low-frequency 
electromagnetic waves could be emitted. If a non- negligible magnetic field is presented in 
the upstream flow, the synchrotron maser instability develops at the shock front (Langdon, 
Arons & Max 1988). 

One-dimensional particle-in-cell (PIC) simulations of relativistic magnetized shocks in 
electron-positron and electron-positron-ion plasmas (Gallant et al. 1992; Hoshino et al. 1992) 
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demonstrate intensive electromagnetic waves both upstream (precursor) and downstream of 
the shock; these waves are generated by synchrotron maser instability at the shock front. It 
was shown that the precursor takes a few-per cent of the flow energy. One can anticipate 
that the upstream flow could be significantly affected by interaction with the precursor; this 
in turn could significantly influence the structure of the shock. For example, absorption of 
the precursor would result in strong deceleration and heating of the flow even though the 
energy of the precursor is small compared with the total energy of the flow; this follows 
immediately from the conservation of energy and momentum in the ultrarelativistic case. 
Of course true absorption is negligibly weak in the case of interest; the precursor radiation 
should be eventually absorbed via non-linear plasma processes, e.g., induced scattering. In 
this case one can anticipate formation of non-thermal particle distribution already in the 
upstream flow, which could significantly affect particle acceleration process at the shock. 

Investigation of complicated nonlinear absorption of the precursor wave is beyond the 
scope of the present research; this would require multidimensional simulations at a very 
large scale. The aim of the present study is to consider a simpler mechanism of interaction 
between the electromagnetic precursor and the upstream flow. This mechanism does not 
assume absorption; however, it operates only in electron-ion flows. The basic idea is the 
following. Because electrons interact with the waves whereas ions do not, some velocity 
difference arises between the electron and ion flows illuminated by powerful electromagnetic 
waves. This is not associated with absorption; just because electrons experience relativistic 
oscillations in the field of the strong wave, the velocity of their guiding centers decreases. 
As ions proceed with the initial velocity, a difference in the bulk velocities of electron and 
ion fluid arises. Therefore a longitudinal electric field is generated so that the electrons 
are accelerated whereas the ions are decelerated. It will be shown below that in the highly 
relativistic case, the energy equipartition is achieved between the electrons and ions before 
the flow arrives at the shock front. 

The article is organized as follows. In sect. 2, the basic characteristics of the electron 
motion in a strong electromagnetic wave are outlined. The behavior of a homogeneous 
electron-ion flow in a strong electromagnetic wave is studied in sect. 3. In sect. 4, the 
synchrotron maser instability at the shock front is considered. PIC simulations of relativistic 
shocks are presented in sect. 5. The results are discussed in sect. 6. In the Appendix, a 
growth rate of the synchrotron instability of a relativistic, narrow ring is estimated. 
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2. Electron in a strong electromagnetic wave 

Let us first briefly outline motion of an electron in a strong wave. Let the electron move 
along the x axis in positive direction and the wave propagate in the opposite direction. Let 
the wave be monochromatic and polarized in the ^/-direction, E y = E smu>(x + 1), E z — 0. 
The electron has two integrals of motion (e.g. Landau & Lifshitz 1987, pp 113-115, 119). 
Invariance with respect to a shift in the y-direction implies conservation of the y-component 
of the generalized momentum, p y — eA y = const, where p = m e u is the particle momentum, 
A the vector potential of the wave, u the electron 4-velocity, m e the electron mass and e 
the positive quantity equal, in absolute value, to the electron charge. The speed of light is 
taken to be unity throughout this paper. Only the y-component of the vector potential is 
non-zero in the linearly polarized wave, 



A y = (Eq/uj) cos uj(x + t); (1) 



therefore one can write 



where 



u y = a cosu(x + t); (2) 

a = 3 

m e uj 

is the strength parameter of the wave. If a <C 1, the wave is linear, at a > 1 oscillations of 
the electron in the field of the wave become relativistic. 

The second integral of motion, which follows from the equations of motions, is 7 + u x = 
const = 7 (J + v x0 ), where 7 = + u 2 x + u 2 y is the electron Lorentz factor. This yields, 
together with Eq.(2), 

a 2 cos 2 (a; + t) 
2(1 + v )7o 

In what follows, we assume 

l«««7o (5) 

so that the wave is strong (oscillations in the guiding center frame are relativistic) but the 
electron energy does not vary much in the course of oscillation. Making use of Eqs.(2, 4) one 
can easily find the velocity of the electron guiding center v gc = (u x /~f), where the angular 
brackets mean averaging over the wave period. As all the velocities are close to the speed of 
light, one can conveniently use the Lorentz factor of the guiding center frame: 



v^o 



7 9 c = (l-<)-^ = ^. (6) 

One can see that in the strong wave, the velocity of the electron guiding center decreases. 
Consequently, if an electron-ion flow is illuminated by a strong electromagnetic wave, the 
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electrons lag behind the ions. In this case a longitudinal electric field arises and electrons 
will be accelerated whereas ions will be decelerated. 



3. Energy exchange between electrons and ions in a strong electromagnetic 

wave. 

To gain an impression of how the relativistic plasma flow interacts with a strong wave, 
let us consider evolution of a spatially homogeneous relativistic plasma flow in a wave of 
a constant amplitude satisfying the condition (5). Let the wave be linearly polarized and 
propagate towards the plasma flow. 



3.1. Non-magnetized flow. 

The electron equations of motion can be written as 

du v ( dA y dA v \ du x / 8A V \ 

Here A y is the vector potential of the wave (1), E\\ the longitudinal (along the x axis) electric 
field. The longitudinal electric field arises because of mismatch in velocities of electron and 
ions; it satisfies the equation 

-jJ- + 4nNe(V x -v x ) = 0, (8) 

where N is the particle density V the ion velocity. The last is found from the equations of 
motion of ions, which are written similarly to Eqs.(7) 

dU y _ f dA y dA y \ _ (IV x _ ( 8Ay 



Here rrii is the ion mass, U = Vr the ion 4-velocity, T the ion Lorentz factor. 

It follows immediately from the first Eq.(7) that Eq.(2) remains valid in this case. 
Combining Eqs.(7) yields the equation 

m e7 ^ = -e£?||5; (10) 

where S = (u x + j)/2. Thus u x + 7 is not constant in the presence of a longitudinal field 
however it varies slowly as compared with the wave period and therefore one can average 
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Eq.(lO) over the wave period to get 

m e ^ = -e(E ll ). (11) 

Here it is taken into account that the right inequality in the condition (5) leads to the relation 
S w (7) (see Eq.(4). 

Similarly the ion equations of motion (9) are reduced to 

U y = —a — - cosuj(x + t); (12) 

Uli 

m *^r = e(i?||); (13) 

which are equivalent to Eqs.(2) and (11). 

Averaging Eq.(8), one can see that (En) is determined by the velocities of the particle 
guiding centers. Making use of the identity 7 2 = 1 + u\ + u 2 , which is conveniently written 
as 25(7 — u x ) = 1 + u^, and Eq.(2), one can write (cp. Eq.(6)) 

<<"> = 1 " < W> = 1 " (14) 

It was taken into account that under the condition (5), S is close to 7 and both are close to 
(7). Analogously one gets for the velocity of the ion guiding center 

/T/\ 1 Z + tamJmj) 2 
Now one can write Eqs.(8), (11) and (13) in the dimensionless form: 

£ + ^- i±i iy^ 0i (16) 
<M +a ^, (17) 

^l-^ w = 0. (18) 

Here r = u p t, uo p = (Ane 2 N /m e ) l/2 , w = (E||)(47rnm e )~ 1/2 . It follows from Eqs. (17) and 
(18) that TOj(r) + m e (7) = (m; + m e )T , where T is the initial flow Lorentz factor. Then 
eliminating w from Eqs. (16) and (17) yields 

d 2 (l) Q 2 2m 2 + a 2 m 2 e = 

dr 2 4(7)2 + 4[(m i + m e )r -m e ( 7 )]2 " 1 ) 
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This equation describes nonlinear oscillations of the electron-ion plasma. The equilibrium 
occurs at 

a/2 + (am e /mi) 2 + am e /mi \ (m;/2m e )r ; a^>mi/m e . 

when the electron and ion guiding centers move with the same velocity (see Eq.(6)). Writing 
down the first integral 

d(l)\ 2 ^ fl2 _|_ 2m 2 + a 2 ml 2m i + 'm e a 2 

dr J 2(7) 2m e [(mi + m e )T - m e {i)\ 2m e T 

one can see that the electron Lorentz factor oscillates around this equilibrium from the initial 
value T up to min(a 2 , M/m)T . So if a > a/ M/m, the energy is efficiently transferred from 
ions to electrons. The characteristic time of the energy exchange, T, is estimated as the 
reverse frequency of oscillations 

I = { a_1/2p o 3 %; a<mi/m e ; ^ 

T j a(m e /miT )~ 3 / 2 Ld p ; a ^> mi/m e . 

If the wave is switched on slowly, so that a grows on a time-scale large compared with T, 
the electron remains in equilibrium and moves with the Lorentz-factor of Eq.(20). Then the 
electron energy grows linearly with a until equilibrium with ions is achieved at a ~ mi/m e . 



3.2. Magnetized flow 

The mechanism outlined works also in magnetized flows. Note that even if the magnetic 
field is tangled in the proper plasma frame, the transverse magnetic field dominates in the 
lab frame therefore it will suffice to consider the flow with the magnetic field B z. Initially 
the magnetic field is frozen into the plasma therefore the electric field v B y is presented 
in the lab frame. In the spatially homogeneous case, the magnetic field remains constant; 
however, the electric field varies according to 

dE . . 

— = Airenvy. (23) 

The equations of motion are obtained by adding the Lorentz force 

F = -e 7 (wy + v x z)B ; (24) 

to the right hand side of Eqs. (7) and the corresponding force to Eqs.(9); here u = E/B is 
the drift velocity. Taking into account that the wave frequency is high, one can introduce the 
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slowly varying variables S = (7 + u x )/2 and R = u y — a cosu(x + 1), which become constant 
in the absence of external fields. Substituting these variables into the electron equations of 
motion and averaging yields 

^f^l±^- 1+u , (25) 

^L^L = R + S £- (26) 
ujb dr 

where ujb = eB /m, £ = E\\/B . Substituting the electric field E = uB into Eq.(23) and 
averaging yields 

ujs_du = _R ^ 
uj p dr S 

Equations (25), (26) and (27) govern the evolution of the electron flow in the presence of the 
magnetic field. 

Let us first neglect the electron-ion coupling and put £ = 0. Then Eqs.(26) and (27) 
yield 

r = 1+2 K) 2 ^ < 28 > 

where r = 21^(1 — u), s = S/Tq. One can check a posteriori that R a. Then Eq.(25) is 
reduced, with account of Eqs.(26) and (28), to 

roS^f^-iVm, (29) 



'dr 2 2T 2 cuj \2s 2 

This equation describes nonlinear oscillations of electrons caused by discrepancy between the 
velocity of electron guiding centers and the drift velocity E/B . If the flow magnetization is 
not too high, 

C^«i; (30) 

the electron energy remains nearly constant, s — 1 1, and Eq.(29) reduces to a linear 
equation 

ro ^^ + s_1 -T' (31) 

which describes oscillations with the frequency uj v I y/To around the equilibrium value s = 
1 + C 2 /4- It follows from Eq.(28) that the equilibrium drift velocity is u = 1 — a 2 /(4rQ), 
which coinsides with the velocity of the electron guiding center (see Eq.(6)). Therefore if the 
wave is switched on slowly such that a grows on a time-scale larger than y^To/up, the system 
evolves remaining in the equilibrium s = sq, u = uq so that the magnetic field remains frozen 
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in the electron fluid in the sense that E/B is equal to the velocity of the guiding-center 
frame. In this case the force (24) remains, on average, close to zero and the evolution of 
the system may be described by non-magnetized equations (16), (17) and (18). So at the 
condition (30), the energy exchange between electrons and positrons proceeds as in the non- 
magnetized case. Note that the characteristic time of the electron-ion energy exchange (22) 
significantly exceeds the period of oscillations described by Eq.(31); this justifies neglect of 
the energy exchange in the above analysis of the interaction of the electron flow with the 
magnetic field. 

If the condition (30) is violated, Eq.(29) describes non-linear oscillations around the 
equilibrium at which the zero-electric-field frame coincides with the velocity of the guiding- 
center frame. The equilibrium electron energy grows with ( and eventually reaches maT 
so that the velocity of the zero electric field frame is equal to the initial velocity. Thus the 
electron energy grows in any case until the velocity of the guiding center frame approaches 
the initial velocity. 



4. Synchrotron maser instability at the shock front 

Relativistic shocks are conveniently normalized by the magnetization parameters 

° s = 4^f ; (32) 

where B is the upstream magnetic field, N the upstream plasma number density (both are 
measured in the shock frame), Y the upstream Lorentz factor. The index s refers to the 
plasma species (i for ions, e for electrons). We consider here the case <7j <C 1 when the 
shock is strong; the electron magnetization may be arbitrary. Note that relativistic shocks 
are generally perpendicular in the sense that the magnetic field is directed predominantly 
perpendicular to the shock normal. This is true unless the field is aligned with the shock 
normal to within 1/T in the upstream frame. 

It follows from the above consideration that the energy exchange between the ions and 
electrons could occur upstream of the relativistic shock provided some electromagnetic in- 
stability at the shock front generates strong enough electromagnetic waves propagating in 
the forward direction. Langdon et al. (1988) and Gallant et al. (1992) demonstrated that 
in relativistic electron-positron shocks, synchrotron maser instability generates strong semi- 
coherent electromagnetic waves both upstream and downstream of the shock. Simulations 
of electron-positron-ion plasma also revealed a strong electromagnetic precursor generated 
by electrons and positrons at the shock front (Hoshino et al. 1992). Heavy ions cannot emit 
at the frequency exceeding the plasma frequency. In electron-positron-ion shocks, the ion 
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synchrotron instability develops because low-frequency magnetosonic waves could propagate 
in electron-positron plasma (Hoshino et al. 1992) however these waves do not propagate 
upstream of the shock because their group velocity is less than the flow velocity therefore 
even in this case the electromagnetic precursor is generated only by electrons and positrons. 
In electron-ion plasma, only the electron synchrotron instability is possible. 

Within the shock structure, the directional motion of the particles is converted into 
the rotational motion in the enhanced magnetic field; one can expect that a ring-like struc- 
ture appears in the particle phase space. This is definitely true for ions however the ion 
synchrotron instability is suppressed by electrons. The synchrotron maser instability is pos- 
sible if electrons form a ring or a horseshoe in the phase space. It is not evident that such 
distributions arise within the shock structure because the electron Larmor radius is much 
smaller than the shock width so that motion of electrons is rather complicated. Simulations 
of non-relativistic shocks reveal holes in the electron phase space (Shimada & Hoshino 2000; 
McClements et al. 2001; Hoshino & Shimada 2002; Schmits, Chapman k Dandy 2002a,b); 
these were attributed to localized electrostatic structures formed at the non-linear stage of 
the Buneman instability. Note that earlier Tokar et al. (1986) found in their simulations 
of the non-relativistic electron-ion shock that an extraordinary mode noise was propagating 
away from the shock; Bingham et al. (2003) interpreted this noise as the electron cyclotron 
maser emission from the shock front. 

In relativistic flows, development of the electrostatic instabilities is suppressed neverthe- 
less one can expect that immediately after the electrons enter the shock front, a significant 
fraction of their kinetic energy is converted into the energy of the Larmor rotation in the 
enhanced magnetic field because in the proper electron frame of reference, the magnetic field 
increases on a time-scale less than the proper electron gyroperiod. The electron with the 
Lorentz factor V "sees" the magnetic field B' = sj B 2 — E\ = B/T, where E± = vB is the 
electric field perpendicular to the flow velocity. In the one-dimensional case, E± remains 
constant and therefore B' varies significantly when B varies only by a fraction 1/T 2 . In the 
shock frame, B varies at the scale of the ion Larmor radius, Ri = miT/eB, so the electron 
"sees" strong variation of the magnetic field when it enters the shock by a depth of only 
d ~ Ri/T 2 . In the proper electron frame, the field increases for the time t' = d/^y ~ -Rj/r 3 
therefore if T > (mi/rrie) 1 ^, the electron motion becomes non-adiabatic when it just enters 
the shock. In this case one can expect formation of a ring-like distribution. In the next 
section, PIC simulations will be presented confirming this conjecture. 

Maser emission from the relativistic ring occurs at the harmonics of the relativistic gy- 
rofrequency Qb = ojb/T- The growth rate of the instability is generally a rather complicated 
function of the plasma parameters. In the Appendix, the instability of a highly relativistic, 
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narrow ring 1 < 5T <C T, is considered in two limiting cases, when the relativistic plasma 
frequency fl p = uo p / y/T is much less and much larger than £l B . Note that a e = fl^/fl 2 ; there- 
fore one can write the corresponding limits as o e ^> 1 and a e <^ 1. The maximal growth 
rate and the corresponding frequency of the emitted waves are written in these limits as 

o.i(r/<$rV/ 4 fl B ; " ~~ \ O^" 1 / 4 ^; a e < 1. 1 J 

One sees that unless <r e is too large or too small, the maser instability develops on a time-scale 
of a few Larmor periods and generates emission at frequencies of a few plasma frequencies. 
One can expect that this remains true also in the intermediate regime a e ~ 1. 

It will be shown in the next section that the electron phase space distribution at the 
front of the highly relativistic shock has a ringlike shape and that such shocks do emit strong 
low frequency radiation. The strength parameter of the emitted waves may be estimated 
as follows. The emission frequency was shown to be something larger than the plasma 
frequency therefore one can write uj = i]Q p , where 77 is a factor of the order of a few in a 
very wide range of a e . The amplitude of the wave is expressed via the wave power, which 
may be parametrized as a fraction £ < 1 of the electron upstream energy. Then the strength 
parameter (3) can be expressed as 

a = y^T . (34) 
V 

It follows from this estimate that the precursor from a highly relativistic shock may be 
considered as a strong wave even if the emitted fraction of the electron energy is small. 
Therefore one can anticipate, according to the results of sect. 3, an efficient energy exchange 
between the electrons and ions upstream of the shock. 



5. Energy exchange between electrons and ions upstream of the relativistic 

shock. 

In order to check whether the mechanism outlined is operative at relativistic shocks, 
we performed PIC simulations. We used a one-dimensional, relativistic, electromagnetic 
code essentially the same as described by Birdsall & Langdon (1991). The simulation is 
one-dimensional in space along the x axis, which coincides with the flow direction. Particle 
motion is restricted to x — y plane, the magnetic field is oriented in the direction z. The 
particles are advanced in time using the relativistic Lorentz force equation. The particle 
momentums, not velocities, were used as variables in order to avoid losing accuracy when 
calculating expressions such as (1 — v 2 )~ l l 2 . The longitudinal electric field is found from 
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the Poisson equation. The transverse fields are evolved via Maxwell's equations reduced to 
equations for the right-hand and left-hand propagating waves, F ± = E±B. These equations 
are solved along the exact vacuum characteristics, x ± t, such that vacuum waves propagate 
exactly with the speed of light; then even a highly relativistic plasma flow does not generates 
numerical Cerenkov emission (Birdsall & Langdon 1991). 

At t — the computation box is filled with a homogeneous plasma moving to the right 
with the Lorentz factor T = (1 — Vq )^ 1 ^ 2 . The thermal velocity dispersion is negligible. The 
plasma initially carries a uniform magnetic field, B Q , as well as the electric field E x = VqBq. 
The flow moves to the right against a wall from which the particles are elastically reflected. 
The plasma is continuously injected from the left boundary. The condition of no incoming 
waves was adopted at the right boundary. In Figs. 1 and 2, the simulation results are 
presented for a run with the initial Lorentz factor of the flow T = 50 and the ion-to-electron 
mass ratio mj/m e = 200. The flow is weakly magnetized, <7j = 0.003. The system grid size 
is 204, 000, the particle density is 2 per species per cell, the time-step, At = Ax, was chosen 
such that eB At/mc = 0.5. The length is measured in units of the upstream relativistic ion 
Larmor radius, Ri = MT Q /eB Q . The physical system size is 8.6-Rj. 

According to the simulation results, a shock arises initially at the right boundary and 
then moves to the left. Fig. 1 shows the initial stage when the shock was just formed. 
One can see a well-developed ion loop in Fig. la. The front of the shock occurs when the 
upstream flow meets the first reflected ions (x = 8. 01 Ri at the presented snapshot). The 
electron flow near this point is shown in Fig. lb; one can see that the electrons enter the 
shock non-adiabatically and that their motion becomes rotational. The ring formed in the 
electron phase space is clearly seen in Fig. lc. Just below this initial region, the phase 
space hole disappears because the synchrotron instability develops and electron distribution 
spreads. Waves generated by the instability are clearly seen in Fig. Id. The electromagnetic 
precursor propagates ahead of the shock front exciting oscillations of electrons and making 
the electrons lag behind the ions. A large-scale longitudinal electric field arises, which 
accelerates electrons. This acceleration is seen in Fig. le. The electrons start to accelerate 
when the precursor reach them (at x = 6.8i?j in the presented snapshot). At this stage, the 
electron energy gain is not large and the electrons enter the shock with an energy still much 
less than the energy of ions. 

In the course of time, the shock propagates farther to the left and so does the precur- 
sor. As the electrons with a larger energy enter the shock, the amplitude of the precursor 
increases. In the wave with a larger amplitude, the electron guiding center slows down, a 
larger longitudinal electric field develops and, according to Eq.(20), the electrons acquire 
larger energy. When these new electrons enter the shock, they emit even stronger waves 
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and then the energy of electrons in the upstream flow increases even more. The ion energy 
decreases appropriately because the total energy of the flow is conserved. Eventually, long- 
wavelength electrostatic oscillations develop. In these oscillations, the average energies of 
electrons and ions are equal. Fig.2 shows the simulation results at t — 8.2i?j. Here the shock 
front is at x = 5i?j. The shock is sharp, the shock width is of the order of Ri. Strong high 
frequency oscillations of the magnetic field at 0.5 < z < 5 (see Fig. Id) represent the electro- 
magnetic precursor. Downstream of the shock, the plasma density and the mean magnetic 
field are 3 times those upstream as they should be in a two-dimensional relativistic gas with 
a ratio of specific heats 3/2 (note that all quantities are measured in the downstream frame; 
in the shock frame, the density jump is 2). 

These simulations show that the energy equipartition between electrons and ions is 
achieved before the plasma enters the shock front so that downstream of the shock, the 
temperatures of ions and electrons are equal. The front of the precursor was emitted when 
the electrons with the initial Lorentz factor entered the shock; therefore the wave amplitude is 
initially not very large and the electrons are accelerated slowly but the amplitude grows when 
the electrons with larger energies enter the shock and therefore the precursor amplitude grows 
and the electrons are accelerated more rapidly. When the electrons achieve equipartition with 
ions, the amplitude of the precursor becomes maximal. At this stage, the precursor is taking 
few per-cent of the total energy of the flow. 

The electrons at the shock front emit electromagnetic waves both forward and backward. 
In weakly magnetized flows, (Jj < 0.01, the energy density of the waves downstream of the 
shock exceeds that of the shock-compressed background magnetic field. 

We ran simulations with the Lorentz factor of the flow from 2 to 50. It was found that 
electrons efficiently take energy from ions provided the flow is highly relativistic, T > 10. 
In the mildly or non relativistic case, the shock does not emit waves strong enough to make 
electrons lag behind ions, so no longitudinal electric field arises. 

There is no sign of nonthermal particle acceleration in these simulations; downstream 
of the shock, the distribution functions of both electrons and ions are Maxwellian. One of 
the reasons is that in the one-dimensional simulations, the electric field is orthogonal to the 
magnetic field. The Fermi acceleration presumably occurs at much larger scale and anyway 
requires three-dimensional turbulence. However this study shows that the nonthermal tail 
in the electron spectrum would begin, provided it forms, from the energy of the order of the 
upstream ion energy. 
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6. Discussion 

It has been demonstrated in this paper that the maser instability develops at the front of 
the ultrarelativistic shock and generates electromagnetic waves propagating both upstream 
and downstream of the shock. Interaction of the upstream plasma flow with these waves 
results in efficient energy exchange between ions and electrons so that the energy equilibration 
occurs in the upstream flow before it arrives at the shock front. An important point is that 
the energy of the maser radiation is generally proportional to the energy of electrons entering 
the shock. Therefore eventually a few percent of the total energy of the flow is radiated away 
in the form of low-frequency waves. 

Non-linear interactions of these waves with the plasma particles could significantly af- 
fect the process of particle acceleration. This issue demands multi-dimensional study and 
therefore is beyond the scope of the present work. In any case, electron-ion equilibration 
facilitates any acceleration process. 

If the magnetization of the upstream flow is low, a < 10~ 2 , the energy density of waves 
downstream of the shock exceeds that of the shock-compressed background field. Relativistic 
particles could radiate in the field of the waves via nonlinear Compton scattering (e.g., 
Melrose 1980, pp. 136-141). The power and characteristic frequencies of this emission 
are similar to those for synchrotron emission in the magnetic field of the same strength 
as the wave amplitude. Therefore the overall spectrum of radiation from electrons with a 
power-law energy distribution generally mimics the synchrotron spectrum. However, there 
are marked differences. For example in the low-frequency range, uj <C (eEo/mec)^ 2 , the 
radiation spectrum of the nonlinear Compton scattering exhibits a larger slope than the 
synchrotron spectrum, oc uo (He et al. 2003) instead of the customary oc oo 1 ^. 

Note that a good fraction of the gamma-ray bursts (GRBs) exhibit X-ray spectra grow- 
ing faster than uo 1 ^ (Preece et al. 1998). There is still no conventional explanation of this 
fact (Lloyd k Petrosian 2000; Medvedev 2000, 2006; Meszaros & Rees 2000; Fleishman 2006) 
The presented model could naturally account for this behavior if one attributes the GRB 
emission, at least partially, to the nonlinear Compton scattering of low-frequency waves gen- 
erated by the synchrotron maser at the shock front. Note also that observations of GRB 
afterglows suggest (Waxman 1997; Frail, Waxman & Kulkarni 2000; Eichler & Waxman 
2005) that the electrons gain a significant fraction of the total energy such that their average 
energy is comparable with that of the ions. The present study shows that electrons indeed 
take about one-half of the total energy, which makes it relatively easy for any acceleration 
process to transfer a significant fraction of this energy to a nonthermal tail. 

I am grateful to David Eichler and Anatoly Spitkovsky for valuable discussions. This 
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Appendix. Synchrotron instability 

Let us assume that all electrons rotate in the magnetic field with the Lorentz factor 
r ^> 1. Such a plasma is evidently unstable with respect to maser synchrotron emission. 
The growth rate of the instability is maximal for the electromagnetic wave propagating 
perpendicularly to the magnetic field and polarized perpendicularly to the magnetic field; 
only this wave is considered here. 



6.i. n B ^>n p 

In this case one can neglect the influence of the plasma on the dispersion properties of 
the emitted waves. Then the growth rate of the synchrotron instability at the j-th harmonics 
is easily found as (e.g. Alexandrov, Bogdankevich & Ruhadze 1984, sect. 6. 4) 



n 



2 „-2 3 



1/3 



where Jj is the Bessel function, V = 1 — l/2r 2 the electron velocity. The growth rate is 
maximal at the first harmonics, uo = Vt B = £l p crl^ 2 , and may be written as 



a e > 1. 



(Al) 



6.2. n B <t:n p 

In this case only high harmonics of the gyrofrequency can propagate, and therefore be 
emitted, within the plasma. Then the maser instability develops only in the Rasin-Tsytovich 
range unless the electron ring is unrealistically narrow (McCray 1966; Zheleznyakov 1967). 
Sazonov (1970) and Sagiv & Waxman (2002) considered the maser instability of relativistic 
electrons with isotropic distribution in the case where the dispersive properties of the medium 
are determined by these electrons. In this section, a similar problem is solved for a ring-like 
distribution. 

Making use of the Einstein coefficient method (e.g. Ginzburg 1989, ch.10) one can write 
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the growth rate of the instability as the absorption coefficient with the opposite sign: 

8tt 3 c 2 



K = 



(Ni — Nj) . 



Here Ni is the population of the i-th Landau level. Introducing the distribution function, 
f(p±,P\\), normalized by the condition f fdp = 1 and making use of the energy conservation 
in the emission process, one can write 



P± 



where e is the electron energy, and p±, p» are the components of the momentum perpendic- 
ular and parallel to the magnetic field. The Einstein coefficients, Aj, are straightforwardly 
expressed via the conventional formula for the power of the synchrotron radiation from a 
single electron, q; this yields 



K 



4n 2 cN f df(p±,0) 



or 



p 2 L dp ± ; 



31/2.2 



2c[l + 7 2 (l-n 2 )] 1 /2 u t 
3 



dp± 

roc / ^ 

/ K 5/3 (z)dz + K 2/3 — 

Ju>/u> c \ U c 



(A2) 



uj Bl '[l + 1 \\-n l 



-3/2. 



where Kj is the Macdonald function, n the refraction index of the medium and it was assumed 
that 7 ^> 1 so that the electron radiates exactly in the direction of its motion. Integrating 
equation (A2) by parts, one gets 



K 



V3i 



10 



f(p±,0) 



_d_ 

dpi 



w/uJc 



K 5/3 (z)dz + K2/3 



dpi 



(A3) 



The refraction index of the wave propagating transversely to the magnetic field and 
polarized transversely to the magnetic field is expressed via the transverse dielectric permit- 
tivity as n = yfe, where (e.g. Alexandrov et al. 1984) 



1 + -^ 



Py 



df 



7(0; - kv x ) dp y 

For a narrow relativistic ring with the distribution function 



dp. 



f(P±,P\\) 



47r 2 m3r(<5r)' 



■ exp 



(p ± - m e Yf + p 2 ' 
2{Sff 



{AA) 
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where ST <C T, one finds 



n 2 I k 2 v 2 



k 2 [ V uj 2 w r 2 Vw 2 - k 2 V 2 
In the case w > fi p so that n - 1 C 1, the refraction index takes the usual form 

n 2 

n =1 

In the Rasin-Tsytovich regime, Y 2 {n 2 — 1) ^> 1, equation (A3) yields for the distribution 
function (A4) 

p 



V31W 



k= F 



where 



P CO 

F(z) = zK 5/3 (z) - / K 5/3 (x)dx - (A/3)K 2/3 (z). 

J z 



The function F(z) achieves its maximum at z — 2.2 and its maximal value is 0.085. Now 
one can estimate the maximal growth rate of the instability as 



The corresponding emission frequency is u = 0.5Q p a t 



-1/4 
e 
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Fig. 1. — Electron and ion phase space and electromagnetic fields at t — 1.9Rf, the initial 
Lorentz factor T = 50, the magnetization cr, = 0.003, the ion-electron mass ratio M/m = 
200. a) longitudinal momentum of ions; b) longitudinal momentum of electrons at the front 
of the shock; c) electron phase space at the front of the shock, 8.088 < x/Ri < 8.101; d) 
magnetic field; e) longitudinal momentum of electrons 
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Fig. 2. — Electron and ion phase space and electromagnetic fields at t — 8.2i?j. Parameters 
are the same as in Fig.l. a) longitudinal momentum of ions; b) longitudinal momentum of 
electrons; c) transverse momentum of electrons; d) magnetic field; e) longitudinal electric 
field 



